Tumor morphological evolution: directed migration and gain and loss of the self-metastatic phenotype

<p>Abstract</p> <p>Background</p> <p>Aside from the stepwise genetic alterations known to underlie cancer cell creation, the microenvironment is known to profoundly influence subsequent tumor development, morphology and metastasis. Invasive cluster formation has been assumed to be dependent on directed migration and a heterogeneous environment - a conclusion derived from complex models of tumor-environment interaction. At the same time, these models have not included the prospect, now supported by a preponderance of evidence, that only a minority of cancer cells may have stem cell capacity. This proves to weigh heavily on the microenvironmental requirements for the display of characteristic tumor growth phenotypes. We show using agent-based modeling that some defining features of tumor growth ascribed to directed migration might also be realized under random migration, and discuss broader implications for cause-and-effect determination in general.</p> <p>Results</p> <p>Considering only the properties of random migration in tumors composed of stem cells and committed cells, we are able to recapitulate a characteristic clustering feature of invasive tumor growth, a property we attribute to "self-metastatic" growth. When the additional influence of directed migrations under chemotactic environments are considered, we find that tumor growth and invasive morphology are supported while the tumor is distant from the source, but are progressively discouraged as the tumor converges about that source.</p> <p>Conclusions</p> <p>We show that invasive clustering can derive from basic kinetic assumptions often neglected in more complex models. While higher-order mechanisms, e.g. directed migration upon chemotactic stimuli, may result in clustering growth morphologies, exclusive attributions of this phenotype to this or other structured microenvironments would be inappropriate, in light of our finding these features are observable in a homogeneous environment. Furthermore, directed migration will result in loss of the invasive phenotype as the tumor approaches the attractor source. Reviewers: This article was reviewed by Mark Little and Glen Webb.</p

Non-stem cancer cell kinetics modulate solid tumor progression

<p>Abstract</p> <p>Background</p> <p>Solid tumors are heterogeneous in composition. Cancer stem cells (CSCs) are believed to drive tumor progression, but the relative frequencies of CSCs versus non-stem cancer cells span wide ranges even within tumors arising from the same tissue type. Tumor growth kinetics and composition can be studied through an agent-based cellular automaton model using minimal sets of biological assumptions and parameters. Herein we describe a pivotal role for the generational life span of non-stem cancer cells in modulating solid tumor progression <it>in silico</it>.</p> <p>Results</p> <p>We demonstrate that although CSCs are necessary for progression, their expansion and consequently tumor growth kinetics are surprisingly modulated by the dynamics of the non-stem cancer cells. Simulations reveal that slight variations in non-stem cancer cell proliferative capacity can result in tumors with distinctly different growth kinetics. Longer generational life spans yield self-inhibited tumors, as the emerging population of non-stem cancer cells spatially impedes expansion of the CSC compartment. Conversely, shorter generational life spans yield persistence-limited tumors, with symmetric division frequency of CSCs determining tumor growth rate. We show that the CSC fraction of a tumor population can vary by multiple orders of magnitude as a function of the generational life span of the non-stem cancer cells.</p> <p>Conclusions</p> <p>Our study suggests that variability in the growth rate and CSC content of solid tumors may be, in part, attributable to the proliferative capacity of the non-stem cancer cell population that arises during asymmetric division of CSCs. In our model, intermediate proliferative capacities give rise to the fastest-growing tumors, resulting in self-metastatic expansion driven by a balance between symmetric CSC division and expansion of the non-stem cancer population. Our results highlight the importance of non-stem cancer cell dynamics in the CSC hypothesis, and may offer a novel explanation for the large variations in CSC fractions reported <it>in vivo</it>.</p

A proliferation saturation index to predict radiation response and personalize radiotherapy fractionation

BACKGROUND: Although altered protocols that challenge conventional radiation fractionation have been tested in prospective clinical trials, we still have limited understanding of how to select the most appropriate fractionation schedule for individual patients. Currently, the prescription of definitive radiotherapy is based on the primary site and stage, without regard to patient-specific tumor or host factors that may influence outcome. We hypothesize that the proportion of radiosensitive proliferating cells is dependent on the saturation of the tumor carrying capacity. This may serve as a prognostic factor for personalized radiotherapy (RT) fractionation. METHODS: We introduce a proliferation saturation index (PSI), which is defined as the ratio of tumor volume to the host-influenced tumor carrying capacity. Carrying capacity is as a conceptual measure of the maximum volume that can be supported by the current tumor environment including oxygen and nutrient availability, immune surveillance and acidity. PSI is estimated from two temporally separated routine pre-radiotherapy computed tomography scans and a deterministic logistic tumor growth model. We introduce the patient-specific pre-treatment PSI into a model of tumor growth and radiotherapy response, and fit the model to retrospective data of four non-small cell lung cancer patients treated exclusively with standard fractionation. We then simulate both a clinical trial hyperfractionation protocol and daily fractionations, with equal biologically effective dose, to compare tumor volume reduction as a function of pretreatment PSI. RESULTS: With tumor doubling time and radiosensitivity assumed constant across patients, a patient-specific pretreatment PSI is sufficient to fit individual patient response data (R(2) = 0.98). PSI varies greatly between patients (coefficient of variation >128 %) and correlates inversely with radiotherapy response. For this study, our simulations suggest that only patients with intermediate PSI (0.45–0.9) are likely to truly benefit from hyperfractionation. For up to 20 % uncertainties in tumor growth rate, radiosensitivity, and noise in radiological data, the absolute estimation error of pretreatment PSI is <10 % for more than 75 % of patients. CONCLUSIONS: Routine radiological images can be used to calculate individual PSI, which may serve as a prognostic factor for radiation response. This provides a new paradigm and rationale to select personalized RT dose-fractionation

A Multicompartment Mathematical Model of Cancer Stem Cell-Driven Tumor Growth Dynamics

Tumors are appreciated to be an intrinsically heterogeneous population of cells with varying proliferation capacities and tumorigenic potentials. As a central tenet of the so-called cancer stem cell hypothesis, most cancer cells have only a limited lifespan and thus cannot initiate or re-initiate tumors. Longevity and clonogenicity are properties unique to the subpopulation of cancer stem cells. To understand the implications of the population structure suggested by this hypothesis - a hierarchy consisting of cancer stem cells and progeny non-stem cancer cells which experience a reduction in their remaining proliferation capacity per division - we set out to develop a mathematical model for the development of the aggregate population. We show that overall tumor progression rate during the exponential growth phase is identical to the growth rate of the cancer stem cell compartment. Tumors with identical stem cell proportions, however, can have different growth rates, dependent on the proliferation kinetics of all participating cell populations. Analysis of the model revealed that the proliferation potential of non-stem cancer cells is likely to be small to reproduce biologic observations. Furthermore, a single compartment of non-stem cancer cell population may adequately represent population growth dynamics only when the compartment proliferation rate is scaled with the generational hierarchy depth

Illuminating the Numbers: Integrating Mathematical Models to Optimize Photomedicine Dosimetry and Combination Therapies

Cancer photomedicine offers unique mechanisms for inducing local tumor damage with the potential to stimulate local and systemic anti-tumor immunity. Optically-active nanomedicine offers these features as well as spatiotemporal control of tumor-focused drug release to realize synergistic combination therapies. Achieving quantitative dosimetry is a major challenge, and dosimetry is fundamental to photomedicine for personalizing and tailoring therapeutic regimens to specific patients and anatomical locations. The challenge of dosimetry is perhaps greater for photomedicine than many standard therapies given the complexity of light delivery and light–tissue interactions as well as the resulting photochemistry responsible for tumor damage and drug-release, in addition to the usual intricacies of therapeutic agent delivery. An emerging multidisciplinary approach in oncology utilizes mathematical and computational models to iteratively and quantitively analyze complex dosimetry, and biological response parameters. These models are parameterized by preclinical and clinical observations and then tested against previously unseen data. Such calibrated and validated models can be deployed to simulate treatment doses, protocols, and combinations that have not yet been experimentally or clinically evaluated and can provide testable optimal treatment outcomes in a practical workflow. Here, we foresee the utility of these computational approaches to guide adaptive therapy, and how mathematical models might be further developed and integrated as a novel methodology to guide precision photomedicine

Cancer Stem Cell Plasticity as Tumor Growth Promoter and Catalyst of Population Collapse

It is increasingly argued that cancer stem cells are not a cellular phenotype but rather a transient state that cells can acquire, either through intrinsic signaling cascades or in response to environmental cues. While cancer stem cell plasticity is generally associated with increased aggressiveness and treatment resistance, we set out to thoroughly investigate the impact of different rates of plasticity on early and late tumor growth dynamics and the response to therapy. We develop an agent-based model of cancer stem cell driven tumor growth, in which plasticity is defined as a spontaneous transition between stem and nonstem cancer cell states. Simulations of the model show that plasticity can substantially increase tumor growth rate and invasion. At high rates of plasticity, however, the cells get exhausted and the tumor will undergo spontaneous remission in the long term. In a series of in silico trials, we show that such remission can be facilitated through radiotherapy. The presented study suggests that stem cell plasticity has rather complex, nonintuitive implications on tumor growth and treatment response. Further theoretical, experimental, and integrated studies are needed to fully decipher cancer stem cell plasticity and how it can be harnessed for novel therapeutic approaches

The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity

The synergy of radiation and the immune system is currently receiving significant attention in oncology as numerous studies have shown that cancer irradiation can induce strong anti-tumor immune responses. It remains unclear, however, what are the best radiation fractionation protocols to maximize the therapeutic benefits of this synergy. Here, we present a novel mathematical model that can be used to predict and dissect the complexity of the immune-mediated response at multiple tumor sites after applying focal irradiation and systemic immunotherapy. We successfully calibrate the proposed framework with published experimental data, in which two tumors were grown in mice at two spatially-separated sites from which only one was irradiated using various radiation fractionation protocols with and without concurrent systemic immunotherapy. The proposed model is calibrated to fit the temporal dynamics of tumor volume at both sites and can predict changes in immune infiltration in the non-irradiated tumors. The model was then used to investigate additional radiation fractionation protocols. Model simulations suggest that the optimal radiation doses per fraction to maximize anti-tumor immunity are between 10 and 13 Gy, at least for the experimental setting used for model calibration. This work provides the framework for evaluating radiation fractionation protocols for radiation-induced immune-mediated systemic anti-tumor responses

Range-Bounded Adaptive Therapy in Metastatic Prostate Cancer

Adaptive therapy with abiraterone acetate (AA), whereby treatment is cycled on and off, has been presented as an alternative to continuous therapy for metastatic castration resistant prostate cancer (mCRPC). It is hypothesized that cycling through treatment allows sensitive cells to competitively suppress resistant cells, thereby increasing the amount of time that treatment is effective. It has been proposed that there exists a subset of patients for whom this competition can be enhanced through slight modifications. Here, we investigate how adaptive AA can be modified to extend time to progression using a simple mathematical model of stem cell, non-stem cell, and prostate-specific antigen (PSA) dynamics. The model is calibrated to longitudinal PSA data from 16 mCRPC patients undergoing adaptive AA in a pilot clinical study at Moffitt Cancer Center. Model parameters are then used to simulate range-bounded adaptive therapy (RBAT) whereby treatment is modulated to maintain PSA levels between pre-determined patient-specific bounds. Model simulations of RBAT are compared to the clinically applied adaptive therapy and show that RBAT can further extend time to progression, while reducing the cumulative dose patients received in 11/16 patients. Simulations also show that the cumulative dose can be reduced by up to 40% under RBAT. Through small modifications to the conventional adaptive therapy design, our study demonstrates that RBAT offers the opportunity to improve patient care, particularly in those patients who do not respond well to conventional adaptive therapy